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We present a new method to obtain topological insulator Bi 2 Se 3 thin films with a centimeter large 
lateral length. To produce amorphous Bi 2 Se 3 thin films we have used a sequential flash-evaporation 
method at room temperature. Transmission electron microscopy has been used to verify that the 
prepared samples are in a pure amorphous state. During annealing the samples transform into the 
rhombohedral Bi 2 Se 3 crystalline strcuture which was confirmed using X-ray diffraction and Raman 
spectroscopy. Resistance measurements of the amorphous films show the expected Mott variable 
range hopping conduction process with a high specific resistance compared to the one obtained in 
the crystalline phase (metallic behavior). We have measured the magnetoresistance (MR) and the 
Hall effect (HE) at different temperatures between 2 K and 275 K. At temperatures T < 50 K and 
fields B < 1 T we observe weak anti-localization in the MR; the Hall measurements confirm the 
n-type character of the samples. All experimental results of our films are in quantitative agreement 
with results from samples prepared using more sophisticated methods. 


I. INTRODUCTION 

Topological insulators (TI) are a new class of materials, 
which were theoretically predicted to exist in different Bi 
based materialsii^. Angle-resolved photoemission spec¬ 
troscopy (ARPES) has been used to verify the topolog¬ 
ical surface states^^— of such materials. TI have, due to 
strong spin-orbit coupling, an insulating bulk energy gap 
and generate conducting topologically protected gapless 
electronic surface states, which are robust against disor¬ 
der and magnetic impurities^. 

These new materials have very interesting properties 
from the basic physics point of view and for future ap¬ 
plications, for example the observation of the quantum 
spin Hall (QSH) effect^ or the realization of Majorana 
fermions for the application in topological quantum com- 
putatioui^i^. The alloy Bi 2 Se 3 was proposed^i as a three- 
dimensional (3D) TI, and later it was experimentally con¬ 
firmed by ARPESiS that this system possesses a single 
Dirac fermion at the surface. 

Diverse investigations of the transport properties in 
Bi 2 Se 3 were done, e.g. resistance measurements as a 
function of the temperature verify the metallic surface 
state&ii, the Hall effect confirms that Bi 2 Se 3 is a n- 
type carrier and the magnetoresistance shows weak anti¬ 
localization (WAL) ^^d^ . The WAL effect is used as evi¬ 
dence that the electrons moving at the surface of the ma¬ 
terial are time reversal symmetry protected. In a TI the 
WAL effect is a result of the strong spin-orbit coupling, 
which, in the presence of defects, puts the backscatter- 
ing at a minimum when a magnetic field or a magnetic 
impurity is not present. Under the absence of localiza¬ 
tion, the surface of a TI remains metallic and is expected 
to exhibit perfect WAL due to helical surface state o^^d'^ . 
Therefore 2D WAL was used to verify the topological 


nature in several casesi^Al, where the temperature de¬ 
pendence of the coherence length is used to prove that 
the WAL originated from the 2D surface states (4 oc 
for 2D WAL and 4 oc for 3D WAL)i^. 

In the last years tremendous efforts were realized in 
order to produce this kind of materials as thin films, 
nanowires, nanocrystals^^ etc., using different techniques 
such as molecular beam epitaxy (MBE) ^°d^ , metal- 
organic chemical vapour depositions^ including mechani¬ 
cal exfoliationSSiSi or chemical vapour transports^. Nev¬ 
ertheless, the preparation of a large area of Bi 2 Se 3 films 
and the possibility to obtain a variety of morphology by 
the combination with electron beam lithography tech¬ 
niques is still a challenge. 

The route of the amorphous-crystalline phase is a pow¬ 
erful method already used to produce high quality thin 
hlms of quasicrystalline AlCuFeSS, AlPdMiiSl and AlP- 
dReSS samples, whose existence in the phase diagram 
does not allow variations of more than 2% of the com¬ 
ponents from the correct concentration. To produce 
the amorphous films we use a method which allows us 
to obtain centimeter-large homogeneous samples. The 
composition of the target material shows a maximum 
aberration of 1% when compared with the source ma¬ 
terial. This technique is called the sequential flash- 
evaporation^S (SEE) and can be easily combined with 
electron beam lithography to obtain samples in the de¬ 
sired shapes, thickness and can also be used to produce 
hetero structures, e.g. for the study of proximity effects. 


II. EXPERIMENTAL DETAILS 

The initial material used for the thin film preparation 
was commercial Bi 2 Se 3 millimetre-large pieces (4-7 mm) 
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Fig. 1. TEM results, a) structure factor S{k) (the inset shows 
the TEM picture), b) pair correlation function g{r). The inset 
shows the reduced atomic distribution function, the linear fit 
at r < 0.3 nm is used to obtain the atomic number density. 


from Alfa Aesar with a purity of 99.999%. After milling 
the pieces to 300-500 p.m in diameter, the crushed pow¬ 
der was introduced into the flash evaporation device. The 
flash evaporation process was done at room temperature 
under a nominal pressure of Ri 5 x 10“^ mbar. The films 
for the transport measurements and X-ray analysis were 
prepared onto Si-substrates (5 x 5 x 0.5 mm^) coated 
with an amorphous insulating 150 nm thick Si 3 N 4 film. 

For transmission electron microscope (TEM) experi¬ 
ments, additional NaCl crystals substrates were placed 
into the chamber for deposition, which, after dilution 
in bi-distilled water, left the films floating on the wa¬ 
ter surface. Thereafter these films were collected onto 
a commercial TEM Cu-grid. Atomic structure analysis 
of the initial material and the prepared thin films was 
measured by means of X-ray diffraction. For TEM mea¬ 
surements we used the transmission electron microscope 
Philips CM 20 FEG, to obtain the structure factor S{K) 
and the pair correlation function g(r), we have used a 
standard method described in^. 

The resistance measurements were done using the con¬ 
ventional four point method; in this case the 5x5 mm^ 
size films were structured by scratching using a micro¬ 
manipulator. The contacts for low temperature measure¬ 
ments were made out of indium and for the high tempera¬ 
ture measurements we used Pd. The ohmic contacts were 
verified by measuring the TV curves. First low tempera¬ 
ture measurements of the amorphous state were done in 
a standard commercial closed cycle refrigerator system 
with a minimal temperature of 25 K. The electrical resis¬ 
tance was measured using a Keithley 6517A electrometer 



Fig. 2. X-ray diffraction of the crystalline and amorphous 
thin films. The obtained lattice parameters are in good agree¬ 
ment with the literature values (shown in brackets). 


due to the high resistance of the amorphous thin films. 

Low-temperature measurements between 2 K and 
275 K were realized in a commercial cryostat (Oxford) 
implemented with a temperature controller (Linear Re¬ 
search LR-700 Bridge) which allows for a temperature 
stabilization of 5 mK over all temperature range. The 
cryostat is also equipped with a ±8 Tesla superconduct¬ 
ing magnet, which allows to measure the MR and Hall- 
effect. 

Low-noise resistance measurements have been per¬ 
formed using the AC technique (Linear Research LR-700 
Bridge) with ppm resolution. In order to determine the 
crystallization temperature we used the DC technique 
(Keithley 2182 with 2001 Nanovoltmeter and Keithley 
6221 current source) on the previously structured sam¬ 
ples. During the measurements a magnetic field was ap¬ 
plied perpendicular to the sample and current, and its 
value was monitored by a Hall sensor installed close to 
the sample. 

The thickness and stoichiometry of the films were 
measured with Rutherford backscattering spectroscopy 
(RBS)^ and particle induced x-ray emission (PIXE)2^ 
using a 2.0 MeV proton beam. There is also an increasing 
interest in the properties of TI materials when magnetic 
impurities are present. Therefore, the amount of mag¬ 
netic impurities was quantified using PIXE^^ measure¬ 
ments as well. This method provides a high sensitivity, 
for the samples studied here the typical minimum detec¬ 
tion limit is ~ 2|j.g/g. The analysis of the RBS and PIXE 
spectra was performed using XRUMP— CeoPIXE iSi. 
The 2.0 MeV protons used have a penetration depth of 
Ri 30 pm in Bi 2 Se 3 and therefore PIXE gives the impurity 
concentration of all the complete sample if its thickness 
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is less than this value, as it is in our case. 

Raman spectra of the crystalline samples were ob¬ 
tained at room temperature and ambient pressure with 
a Dilor XY 800 spectrometer at 514.53 nm wavelength 
(green) and a 2 pm spot diameter. The incident power 
was varied between 0.5 mW to 3 mW in order to check for 
possible sample damage or laser induced heating effects. 
No damage or significant spectral change was observed 
in this range of incident power. 


III. RESULTS AND DISCUSSION 

After the SFE preparation process at room tempera¬ 
ture, the amorphous samples were obtained. The electron 
diffraction pattern of the amorphous sample is shown in 
the inset of Fig. [TJ After analysis we obtained the atomic 
structure factor S{K) and the pair correlation function 
g{r), both results are plotted in Fig. [TJ 

The structure results are comparable to known simi¬ 
lar amorphous semiconductor systems^, but specially to 
the structure of the liquid phase of Bi 2 Se 3 investigated 
by Usuki^. Our results show similar features as the liq¬ 
uid phase, e.g. the position of the first and second peak 
in the S{K) and g(r) are almost at the same values. The 
first peak at ri ~ 0.29 nm in g(r) corresponds to the same 
value in the liquid phase and is close to the corresponding 
first nearest neighbor distance in the crystalline phase of 
Bi 2 Se 3 ^; the peak r 2 ~ 4.18 nm was identified in the liq¬ 
uid phase as consequence of the building of a short-range 
order related to the crystalline phased. After integra¬ 
tion, the first peak in g(r) allows us to calculate the co¬ 
ordination number 3.2. This value is less than that 
of the liquid state (Ac^Liquid ~ 4.12)^ and is due to the 
amorphous phase, which builds more defined first neigh¬ 
bors in comparison to the liquid phase. From the slope of 
the reduced atomic distribution G'(r) = 47rrno[g(r) — 1] 
at small values of r (dashed line in Fig.JTJ), we can calcu¬ 
late the atomic number density no = 35.05 nm“^, which 
is comparable to the value of the liquid phase when cal¬ 
culated using a simple linear extrapolation of the density 
of Bi and Se in the liquid state. From the structure re¬ 
sults and the high temperature stability of the amorphous 
phase (discussed later), we can conclude, according to 
the well established interpretatio n^^i^^ , that the system 
is a stabilized electron-phase, which means the forma¬ 
tion of a pseudogap at the Fermi energy^ is responsi¬ 
ble for the high stability. Consequently, the influence on 
the electronic transport yields a very high resistivity and 
semiconducting-like behavior as a function of the temper¬ 
ature. After annealing in argon atmosphere to 425 K, the 
sample transforms into the crystalline phase. In Fig. |2| 
we present X-ray results of the amorphous and crystal¬ 
lized thin films. In case of the crystalline films we can 
recognize the peaks corresponding to the rhombohedral 
i?3m (T?!^) structure of the Bi 2 Se 3 phase. The calcu¬ 
lated lattice parameters are in very good agreement with 
literature data^^ within experimental precision. Using 



Temperature T(K) 

Fig. 3. Resistance vs. temperature of the amorphous sample 
after preparation at 320 K. The inset shows the fit and the 
two contributions mentioned in the text. 


the Scherrer equation we can determine the crystalline 
grain size, for that we used the most intense peak and 
the calculation gives a value of ss 19 ± 2.2 nm. Taking 
into account that the unit cell is composed of three layers 
stacked along the c-axis, each known as quintuple layer 
(QL), and that a QL has a thickness of ~ 0.96 nm, our 
films contain grains which include 20 QLs. That is im¬ 
portant because in order to have a Dirac cone it is neces¬ 
sary to have at least 5-6 Ql4i. Resistance measurements 
of the as-prepared amorphous films were also done, the 
results are shown in Fig. |3| In our device the highest 
measurable value of the resistance has been reached at 
30 K. In the complete temperature range the resistance 
shows an thermally mechanism, but with two different 
regimes, which can be seen in two different slopes shown 
in the inset of Fig. |3| We have found that the temper¬ 
ature dependence of the resistance can be well described 
by two different variable-range hopping (VHR) conduc¬ 
tion processes. The general equation for VRH is given 
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Temperature T(K) 

Fig. 4. Resistance vs. temperature from 2 K to 275 K of 
the crystalline Bi 2 Se 3 sample. The upper left inset shows the 
resistance when the sample transforms from amorphous to 
crystalline phase, the intersection of the dash lines gives the 
crystallization temperature Tk. The inset on the r.h.s. shows 
the resistance of the crystalline sample. 


by 


R{T) = Rq ■ exp 



( 1 ) 


where Ro is a prefactor and Tq is a characteristic temper¬ 
ature coefficient. The value of the exponent p depends 
on the nature of the hopping process. We found that in 
our amorphous Bi 2 Se 3 the value of p = 0.25 describes 
very well the measured transport behavior. This value 
of p was derived by Motl4^, and that is in the case of a 
constant density of states at the Fermi level. According 
to Mott, Tq is a coefficient which depends on the density 
of states N{Ef) at the Fermi level in the form; 


To = 


18 

kBeNiEp)^ 


( 2 ) 


where ^ is the localization length and ks is the Boltz¬ 
mann constant. As we mentioned before, it is evident 
that the resistance in our sample shows two different 



Applied Field PqH(T) 

Fig. 5. Magnetoresistance vs. applied magnetic field, mea¬ 
sured at different constant temperatures. 


temperature regions, in each one the resistance can be 
described by the Mott-VRH. In order to describe the 
total temperature dependence we used a simple circuit 
containing two resistors in parallel which result is shown 
in Fig. [3]as a line. A deviation of the fit from the experi¬ 
mental data at temperatures T > 320 K can be observed 
because the sample was prepared at ambient temperature 
and above this value the resistance starts to decrease ir¬ 
reversibly. The contributions of the two resistances are 
evident in the inset. From the fit we obtain the character¬ 
istic temperature Tq at every region. For the low temper¬ 
ature region To_lt = 1-02 x 10® K and for the high tem¬ 
perature region Tq-ht = 75.87 x 10® K. Using these val¬ 
ues we can estimate the localization length ^ if the value 
of N{Ef) is known, but this is not yet available in the 
literature. Assuming a value of the same order for sim¬ 
ilar amorphous systems, N{Ef) ~ 10^®^^^)eV“^ cm“®, 
we obtain at low temperatures ^lt ~ 0.5... 3 nm and 
at high temperatures ^ht ~ 0.14.. .0.6 nm. These val¬ 
ues correspond approximately to the first and second 
nearest neighbor distances in the pair correlation func¬ 
tion obtained from the TEM investigations. By means 
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Fig. 6. Low temperature magnetoresistance, the lines show 
the fits using the WAL theory. 



Fig. 7. Coherence length /^(nm) vs. temperature 
y-i/2(K“i/2) Q, (inset) obtained from the fits of the mag¬ 
netoresistance. The red lines in the figures are the fits and 
the expected theoretical values of and a respectively. 


of the VHR analysis it is difficult to conclude on the 
reason of the change in the localization length at spe¬ 
cific temperature ranges which remains an open ques¬ 
tion. The transition from amorphous to crystalline 
phase was monitored by the measurement of the resis¬ 
tance during annealing in a home made device inside an 


oven under Ar-flow. The results are shown in the insets 
of Fig. m The crystallization occurs in a nearly step¬ 
like process and it is evident that the as-prepared film 
has a larger resistance in comparison to the crystalline 
one. The value of the crystallization temperature Tk 
was determined by intersection of the two lines, one be¬ 
fore and another one during the crystallization process. 
We obtain a Tk of « 410 K, which is relatively high in 
comparison to amorphous metallic system o^^d^ but is in 
the same range of amorphous-like materials such as Al- 
CuFe^^i^ and AlPdM n^^d^ , where a pseudogap at the 
Fermi energy was observed by means of angular-resolved 
photoemission studies^ and tunneling spectroscopy^. 

In Fig. |4] we present the results of the temperature 
dependence of the resistance of our film in the crystalline 
phase. At temperatures between 15 K and 285 K the 
sample shows the typical metallic behavior already ob¬ 
served in this system. In the low temperature range 
T < 15 K we observe an increase in the resistance with 
decreasing T, which nearly follows a logarithmic temper¬ 
ature dependence, as observed in some metal o'^^d^ when a 
small amount of magnetic impurities is present. One pos¬ 
sible explanation for this behavior would be an apparent 
Kondo effect^. Note that the resistance increases nearly 
logarithmic below 12K. Atr<3K the slope changes 
which can be a result of the strong influence of the large 
magnetic impurity content to the spin-orbit interaction 
characteristic for this system. 

As consequence of the resistance behavior at low tem¬ 
perature we did some FIXE measurements. We have 
analyzed the initial bulk material and the deposited thin 
film before and after it was transformed from the amor¬ 
phous phase into the crystalline one. The results of the 
elements with the highest concentration are summarized 
in Table HI With the exception of Fe, the concentration 
of all other magnetic elements was below the detection 
limit. According to the FIXE results our preparation 
method has not produced considerable change of the ini¬ 
tial stoichiometry. This result confirms the advantages 
of preparing with the SFE method as before mentioned. 
According to FIXE, the Fe (Cr) content already present 
in the initial material has increased more than one order 
of magnitude after preparation. This contamination can 
only be due to the milling process before SFE, because 
the amount in the amorphous and crystalline are quite 
similar, therefore the contamination in the oven can be 
ruled out. 

In the following we show the results of the transport 
properties under the application of a magnetic field. In 


TABLE 1. Results of FIXE measurements. 


Sample 

Concentration 
(atom %) 

Fe cone. 

(in ppm) 

Cr cone. 

(in ppm) 

Initial mat. 

BUoSeso 

34± 11 

< 29 

Amorphous mat. 

Bi39,4Se60.6 

878 ± 53 

319 ± 58 

Crystallized mat. 

Bi4o,2Ses9.8 

989 ± 48 

339 ± 59 
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Applied Field HqH(T) 


Fig. 8. Hall resistivity vs. applied field at different temper¬ 
atures. The inset on the l.h.s. shows the calculated mobility, 
the inset on the r.h.s. the calculated carrier density. 


spite of the magnetic impurity content and the large poly¬ 
crystallinity of the film, the results are quite similar to 
those reported in the literature for this TI. 

The normalized magnetoresistance (MR) is defined as 


MR = 


R{B) - R{B = 0) 
R{B = 0) ’ 


( 3 ) 


and are plotted in Fig. [S] In general the MR is posi¬ 
tive for all measured temperatures and at high fields and 
high temperatures (B > 1.5 T, T > 100 K) it is domi¬ 
nated by a parabola like dependence, due to the Lorenz 
force contribution to the carriers scattering under influ¬ 
ence of a perpendicular magnetic field^. At tempera¬ 
tures T < 75 K and in the low-field regime {B < 1.5 T) 
the MR increases sharply as the magnetic field increases, 
which is the typical behavior of the weak anti-localization 
(WAL) effect^. When applying a magnetic field, the 
time-reversal symmetry breaks down with increasing field 
strength and the electrons begin to backscatter such that 
the resistance rises. In order to fit our results with 
the theory, we define the magnetoconductance (MC) as 
AGiB) = G{B) - G(0), with G(0) = 1/R(0). The low 
field part {B < 1.5 T) of the MC can be well fitted using 



Fig. 9. Raman spectra of the crystalline BLSea sample. The 
dashed lines indicates the characteristic modes. 


the standard WAL theory developed by Hikami-Larkin- 
Nagaoka, given as: 

AGwal(R) = G{B) - G(0) 



where a is a coefficient whose value was predicted to be 
1/2 in the case of systems with strong spin-orbit interac¬ 
tion, i.e. in the pure WAL regime, ^ is the digamma func¬ 
tion, Bfj, is the dephasing magnetic field related to the 
phase coherence length according to B^j, = h/Ael"^. Us¬ 
ing the Eq. |4l we fit the experimental results, see Fig. |6l 
and estimate the different parameters at each tempera¬ 
ture. The prefactor a and the phase coherence length 
are shown in Fig. [T] According to our fit results, at the 
low temperatures (T < 20 K), a varies around 0.5 ±0.2, 
which is close to the theoretical expected value. The high 
temperature values are between 0.38 and 0.45, this can be 
a consequence of the diminishing of the coherence length 
with rising temperature. The coherence length varies 
from 300 nm (T = 2 K) to 100 nm (T = 50 AT), which 
is in concordance with the value reported for Bi 2 Se 3 thin 
films prepared by molecular beam epitaxy (MBE) with 
a thickness of 20 QL^. According to our results, the 
temperature dependence of the coherence length is T^, 
with b = 0.51 ±0.13. The exponent gives information 
about the dimensionality of the system. According to 
the theory, for a two-dimensional system, 6 = 0.5 and 
for one dimension, b — 1/3. Our result is close to the 
exponent expected for a two-dimensional system^. For 
a and the the observed behavior and values obtained 
are consistent with that reported for thin film o^^i^'^i^^ and 
microwires of Bi 2 Se 3 — . 

The Hall measurements were done at the same tern- 
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peratures as the MR and are presented in Fig. |8l Its 
sign indicates that the major carriers are electrons, i.e. n- 
type over the investigated temperature range. From the 
slope we can calculate the 3D carrier density, the results 
are plotted in the inset of Fig. [8l The value obtained 
n® « 3.7 X 10^®cm“^ is close to that reported in^ 
where the sample was grown without dopants but also 
to^ where Fe was used as dopant. 

According to theoretical predictions, the presence of 
magnetic impurities in a TI should open a small gap in 
the surface states due to the breaking of time-reversal 
symmetry^. To verify these predictions experimentally, 
films containing Fe in unknown concentrations (but with 
less than 2 % according to the authorsiii^) were pro¬ 
duced. For all doped samples the results show the typical 
behavior of the Kondo effect in the temperature depen¬ 
dence of the resistance. The results of the magnetoresis¬ 
tance show the typical WAL behavior and only in samples 
with a high concentration of Fe, a mix of weak localiza¬ 
tion (WL) and WAL effects has been seenii. Anoma¬ 
lous Hall effect was observed in Bi 2 Se 3 thin films when 
doped with Cr, but only when the concentration is > 
1.4 %^. Our Hall effect results do not show any sign 
of an anomalous contribution as a consequence of the 
Fe and Cr content. Similar results were obtained when 
Bi 2 Se 3 was doped with Fe with a concentration of less 
than 2 %^. We conclude that the amount of Fe and 
Cr present in our sample apparently does not affect the 
expected behavior— of the investigated transport prop¬ 
erties of the TI. 

Finally, in the crystalline state we have done Raman 
scattering spectroscopy of our films. According to the 
stoichiometry, the primitive unit cell of Bi 2 Se 3 contains 
five atoms, therefore 15 degrees of freedom are allowed 
at the center of the Brillouin zone {q — 0), from which 
3 are acoustic modes and the rest are optical phonons 
(2Alg, 2Eg, 2 aiu and 2eu)- According to previous worb^l. 
Eg and Aig modes can be distinguished. In our films we 
were able to detect the high frequency Eg and Ajg modes. 


at Rs 133.14 cm“^ and ~ 175.5 cm“^, respectively, see 
Fig. ini The measured Raman shift in our samples is in 
agreement with the results from literature^i^, indicating 
that our films have similar structural quality as others 
produced with more sophisticated methods. 

IV. CONCLUSION 

TI Bi 2 Se 3 polycrystalline thin films can be obtained 
by annealing starting from the amorphous phase. The 
preparation was done using a simple and inexpensive 
preparation method, which can also be used to produce 
similar TI materials, giving the chance to study diverse 
properties before, during and after the disorder-order 
transition. The material in the amorphous phase has a 
high resistivity and behaves like a semiconductor, which 
transforms in to the metallic crystalline phase at temper¬ 
atures higher than room temperature. TEM microscopy, 
X-Ray and RAMAN measurements were used to verify 
that we have produced the Bi 2 Se 3 films in the desired 
crystalline phase. Using FIXE, a method with high ac¬ 
curacy, we have quantified the presence of magnetic im¬ 
purities (such as Fe) in the initial material and in the 
prepared thin films. Nevertheless, despite the films being 
polycrystalline and that magnetic impurities are present, 
all measured electronic transport properties characteriz¬ 
ing a TI like Bi 2 Se 3 were measured and our results are 
similar to those reported in literature. The fact that our 
results are quantitatively comparable with data obtained 
for samples prepared by more sophisticated methods and 
of apparently much better quality, indicates that: either 
the surface states are extraordinarily protected and ro¬ 
bust against disorder and magnetic impurities, or the re¬ 
sults we are comparing are not the ones that one should 
take as characteristics of a good TI. Future systematic 
measurements as a function of disorder and magnetic im¬ 
purities would consequently clarify the amount of disor¬ 
der and magnetic impurities or both, which causes Bi 2 Se 3 
to become as a non- TI material. 
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